Local-spin-density functional for multideterminant density functional theory 
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Based on exact limits and quantum Monte Carlo simulations, we obtain, at any density and spin 
polarization, an accurate estimate for the energy of a modified homogeneous electron gas where elec- 
trons repel each other only with a long-range coulombic tail. This allows us to construct an analytic 
local-spin-density exchange-correlation functional appropriate to new, multideterminantal versions 
of the density functional theory, where quantum chemistry and approximate exchange-correlation 
functionals are combined to optimally describe both long- and short-range electron correlations. 
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I. INTRODUCTION 



Density functional theoiy^'^ (DFT) is by now the 
most popular method for electronic structure calcula- 
tions in condensed matter physics and quantum chem- 
istry, because of its unique combination of low compu- 
tational cost and high accuracy for many molecules and 
solids. There are, however, exceptions to such an ac- 
curacy. Even the best approximations of its key ingre- 
dient, the exchange-correlation (XC) energy functional, 
cannot describe strong electron correlations, like those 
of the cuprates, and cannot exactly cancel the so-called 
self-interaction, a property which the exact functional 
should satisfy. On top of that, they fail to recover 
long-range van der Waals interactions,'* are not com- 
pletely safe for the description of the hydrogen bond;^ 
and have intrinsic problems with situations of near de- 
generacy (when two sets of Kohn-Sham orbitals have very 
close energies) jSi^ more generally, the "chemical accu- 
racy" has not yet been reached. To overcome the lat- 
ter group of problems, there has been a growing inter- 
est in "mixed schemes" which combine the DFT with 
other approximate methods by splitting the coulom- 
bic electron-electron interaction 1/r = Veeif) into a 
short- rang^SR) ^ ^^^^^^it ^sji ^'^ P^^^^ {see^ e.g. 

use different approximations for the LR and the SR con- 
tributions to the exchange and/or correlation energy den- 
sity functionals of the Kohn-Sham (KS) scheme. It de- 
scends from the observation that LR correlations, poorly 
described by local or semi- local density functionals, can 
be accurately dealt with by other techniques, like the 
random-phase approximation (RPA) or standard wave- 
function methods of quantum chemistry. Conversely, 
correlation effects due to the SR part of the electron- 
electron interaction are in general well described by local 
or semilocal functionals>i2iS£ The error function and its 
complement 



erfc(/xr) erf(/^r) 



where fi controls the range of the decomposition, are of- 
ten used^i'^i^i*°i"i*^i*^i*^i^'^ to split the Coulomb interac- 
tion into a SR and a LR part, since they yield analytic 
matrix elements for both Gaussians and plane waves, i.e., 
the most common basis functions in quantum chemistry 
and solid state physics. Correspondingly, the universal 
fimctional of the electron density n, as defined in the 
constrained-search formalismi^i 



F[n] = min(^'|T- 



(2) 



can be divided into a short-range and a complementary 
long-range part F[n] 



mm 



LR 



[n] = F[n] - F^^ln] 



(3) 



or, alternatively, into a long-range and a complementary 
short-range part F [n] = F^r [n] + FgR [n] 



FUn] 



min + KrI*^) 

F[n]-Fi:^[n]. 



(4) 



(1) 



The two decompositions lead to different strategies and 
XC energy functionals, whose merits and drawbacks are 
discussed in Ref. 122- In any event, for actual electronic- 
structure calculations to be performed, these function- 
als ultimately need approximations, in analogy with the 
standard DFT. Regardless of the strategy adopted, the 
potential superiority of "mixed schemes" comes into play 
precisely at this stage: compared to the standard version, 
a DFT which only deals with the SR part of the electron- 
electron interaction should be much more accurately ap- 
proximated, as mentioned, by local-density XC energy 
functionals ji^iiSiSS While both decompositions [Eq. Q 
and Eq. are aimed at the exploitation of the DFT 
scheme for the SR part of the interaction only, the cor- 
responding approximate functionals require an accurate 
description of the homogeneous electron gas (HEG) ei- 
ther with SR [Eq. ©I or LR interaction [Eq. Q]. 
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Up to now, the HEG exchange-correlation energies as 
a function of the cutoff parameter /x and of the elec- 
tron density are available for the SR case from Quan- 
tum Monte Carlo (QMC) simulations, ^'^ and for the 
LR case from coupled-cluster (CC) calculations 7-^^ A 
parametrization of the CC data for the XC energy of the 
HEG with long-range-only interaction has been used in 
Refs. 8 and 16 with very promising results for closed-shell 
systems. Generalized-gradient-corrected density func- 
tional have also been designed and tested within this 
frameworkjiSiii but all existing functionals are limited 
to the spin-unpolarized case. 

The purpose of this paper is to provide, based on novel 
exact limits and quantum Monte Carlo simulations, an 
accurate representation for the energy of the LR-only in- 
teracting HEG not only as a function of the cutoff pa- 
rameter /X and of the total electron density, but also as 
a function of the spin polarization [i.e., as a function 
of the spin densities ^^(r) and ^^(r) separately]. Since 
von Barth and HedinS^ showed, in 1972, that the task 
of finding good approximations to exchange-correlation 
density functionals is greatly simplified if the functional 
is expressed in terms of the spin densities, and that this 
is the simplest way to satisfy the requirement (Hund's 
rule) that a state with larger spin tends to be energet- 
ically favored, the importance of including such a spin 
dependence in approximate functionals was confirmed 
by countless calculations for molecules and solidsi^SiSi In 
this context the decomposition of Eq. |0J, based on the 
constrained-search formalism^Si is generalized to spin- 
DFT as foUowsiSi 



(5) 



The spin-polarized LR-only gas, for which no previous re- 
sults are to our knowledge available, is appropriate, via 
Eq. to a very promising "multideterminantal" ver- 
sion of the spin-DFT. The final outcome of this work 
is thus a local-spin-density approximation of the corre- 
sponding XC functional, given in analytic form, by which 
electronic structure calculations of this new type will be 
possible for unpolarized systemsS^ and, more important, 
for spin-polarized systems, for which no such functional 
is presently available. Such a functional also represents 
the key ingredient for extending gradient-corrected SR 
density functionalsi^iii to spin-DFT. 

The paper is organized as follows. In Sec. [H] we de- 
fine the hamiltonian of the HEG with LR-only interac- 
tion and we derive some exact limits of the correspond- 
ing correlation energy, which is then computed (for val- 
ues of the relevant parameters not accessible to analytic 
methods) with QMC in Sec. Hill The results of Sees. |Il| 
and mil are then used in Sec. llVl to construct an analytic 
parametrization of the LR correlation energy. Sec. 
recalls, calculates and provides in analytic form an alter- 
native definition of the LR correlation energy which in- 
volves the use of pair-correlation functions (also obtained 



from our QMC simulations) and may be of interest within 
optimized-effective-potcntial schemesi^ 

Hartree atomic units are used throughout this work. 



II. DEFINITIONS, AND EXACT LIMITS 

After decomposing the standard (Coulomb interac- 
tion) spin-DFT functional according to Eq. jSJ, the re- 
sulting SR functional ^^sR I'^Ti "-il tie further decom- 
posed, as usually, into a Hartree and an XC term 



EUn]^E^[n] = ^ I dr f dr' n{r) n{r') v^^{\r - r' \) (6) 



Kclnh^l] = ^SR["T''^i] -^hW- (7) 

The local-spin-density (LSD) approximation amounts to 
replacing the exact, unknown functional of Eq. Q with 

■E^xc,LSD[»^T'"i] = (8) 



dvn{r) [e.c(r.(r),C(r)) - eLR.(^,(r), C(r), /i) 



where exc('T-ti"i) is the exchange-correlation energy per 
particle of the standard jellium model^fii^ii^SiSi with uni- 
form spin densities n|,nj, and e^^(n|,n|,^) is the cor- 
responding quantity for a jellium model with LR-only in- 
teraction i'lp^(?'), which forms the object of this paper. In 
the third line of Eq. (jSJ we express the same quantity in 
terms of the electronic density n = n-^+ni = {Airr^ 
and spin polarization ^ = {n]- — ni)/n, thus introducing 
the notation used in what follows. To obtain e^^^r^, C, m) 
we consider a uniform system with LR-only interaction 



N 



-"lr — 



(9) 



i=l 



where V^^i the modified electron-electron interaction 



N 

E 



erf(^|r, -r^l) 



(10) 



is the interaction between the electrons and a rigid, 
positive, uniform background of density n 

N 



eb 



-n 



Y^Jd. 

i=i •' 



erf(^|ri - x|) 



(11) 



and V^^ is the corresponding background-background in- 
teraction 



T/M /"w /"w / erf(/i|x-x^|) 



(12) 
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Our hamiltonian H^^ , and thus its ground-state energy 
per electron e'^^, depends on the density parameter r^, 
on the spin-polarization C,, and on the cutoff parameter /i. 
When /i — > cxD, we recover the standard jellium model; in 
the opposite limit /i — *■ 0, we recover the nonintcracting 
electron gas. In this Section we derive the asymptotic 
behavior for /i — s- and — > oo of the correlation energy 
per electron, defined as 



(13) 



where ts{rsX) — 3fc| tpsiO/^^ is the usual kinetic energy 
of the nonintcracting electron gas, with fcp = (ars)~^, 
a = (4/97r)i/^ and 



1 r 
2 



(i+C)"/' + (i-C)"/' 

the exchange energy is given byi 



(14) 



(15) 



(2y-42/3)e-i/V-3y-H4y3 + 



/tt erf 



(16) 



and has the asymptotic behaviors 



^LR 



27r 



A. Approaching the non-interacting gas 

When /i — > and/or — s- 0, we are approaching the 
limit of the non-interacting Fermi gas. Toulouse et ali^ 
have studied the /z ^ limit of the long-range exchange 
and correlation energy functionals for confined systems 
(atoms, molecules) using standard perturbation theory. 
Their results cannot be applied to the case of an ex- 
tended system like the uniform electron gas, because the 
integrals of their Eqs. (17) and (20) would diverge. In- 
stead, the fi ^ limit (as well as the — > limit) of 
the uniform electron gas can be studied with RPA(2Si24 
which becomes exact both for long-range correlations 
(/i — > 0: in this limit the long-range coulombic tail shows 
up only beyond larger and larger interelectronic distance 
r ~ l//z) and in the high-density hmit (r^ 0). We 
generalize to the LR-only interaction erf (/ir)/r the stan- 
dard RPA expression for the correlation energy (see Ap- 
pendix^for details), and find that, for small [i^Jrl (i.e.. 



small- /z and/or 
scales as 



limit) , the correlation energy e 



LR 



,LR 



= [Mor 



Q{x), 



(19) 



where (f>2{C) is given by Eq. H14() . and the function Q{x) 
has the following asymptotic behaviors 



(20) 



21n(2) 



■ ln(a;) + const. (21) 



Q{x oo) = 



The scaling of Eq. (|19|l for the long-range correlation en- 
ergy was also expected from the fact that the long-range 
part of the pair-correlation function of the standard jel- 
lium model has a similar scaling Notice also that, 
in the small-/i expansion of [Eqs. H19|) and (|20|l ]. the 
term proportional to /i^ exactly cancels with the corre- 
sponding term in the exchange energy e^^ [Eq. (|17|l ]. so 
that the XC energy (exchange plus correlation) has no 
/i^ term; in a confined system, on the other hand, the /i^ 
terms are separately zero for exchange and correlatiouii^ 
We found that the function Q{x) (see Appendix is 
accurately approximated by 

o(.).El^i„( '+;; + '-f + 

TT^ \ I + ax + ax / 

with a = 5.84605, c = 3.91744, d = 3.44851, and b = 
d-37ra/(41n(2)-4). 

A final remark on the scaling of Eq. (|19|l is that, al- 
though it exactly holds only in the small- a; regime, even 
in the large- a; regime of Eq. lpT|l (obtained e.g. with 
small Ts but very large ^), it represents an excellent ap- 
proximation, because in this regime the ^ dependence of 
e^R/Q is described by a function [Eq. (32) of Ref. HI 
which is very similar to [<^2(C)]^j ^^'^ exactly equals it 
for ( — and C = 1- the densities corresponding to 
Vs 0.1 are not affected by this small difference in the C 
dependence of Eq. ((T^ . as discussed in Appendix 1X1 



B. Approaching the coulombic gas 

The large-/! behavior of the long-range correlation 
functional appropriate to the decomposition of Eq. jSJ, 
obtained in Refs.l^andU^ is straigthforwardly extended 
to the uniform electron gas. For C 7^ 1 we hava^S 



,LR 



(rsXilJ-) 



3(0, r„C) 



/27r 11^ 



3gc(0,rs,C) 
8 rf 

o(m-^) 



(23) 



where ec{fs,C) is the correlation energy of the standard 
electron gas with Coulomb interaction, 5(0, r^,^) its on- 
top pair-distribution functionjSSi^LiSi and gc{0,rs,Q = 



4 



g{0,rs,C) — 5(1 — C^)- For the fully polarized gas (C = l) 
the terms proportional to /Lt~^ and in the large-/Lt 
expansion of vanish, and the next leading terms are^ 



ec(r.,C=l)- 



64 r3 ^4 
0(^-«)(24) 



where (7"(0, r^., C= 1) is the second derivative at r = of 
the pair-distribution function'^^-'^^-*" of the fully polarized 
gas, and g'^{0, r„ C= 1) = g"(0, r„ C = 1) - 25/3fc|/5. 



III. DIFFUSION MONTE CARLO 

The details of the simulations are similar to our previ- 
ous calculation of a local density functional for a short- 
range potential. ■^■^ Here we give a technical summary fo- 
cusing on the main differences, which concern the size 
extrapolation and the treatment of the long-range tails 
of the interaction and of the pair pseudopotential in the 
trial wave function. For the reader not keen on techni- 
calities, it is enough to say that we provide a very tight 
upper bound to the exact ground-state energy, choos- 
ing a level of approximation which closely matches the 
Ceperley-Alder^ (GA) result for /i — > 00. 

The ground-state energy of the Hamiltonian of Eq. ^ 
is computed with the diffusion Monte Carlo (DMC) 
method in the fixed-node (FN) approximation,**^ using a 
standard Jastrow-Slater trial function with plane-wave 
orbitals and RPA pseudopotentials42i Several values of 
the density {vg = 1, 2, 5, and 10), of the cutoff parameter 
(/xr^ = 0.5, 1, 2 and 4) and of the spin polarization (C = 
and 1) are considered. The results are fitted (see Sec. lIV|) 
to a convenient analytical expression for the correlation 
energy ec(rs,/i, C)j which also embodies the exact limits 
of Sec. m and is further constrained to recover the CA 
result for the Coulomb potential. 

This constraint sets the target precision of our simu- 
lations, since there is no point in pushing the accuracy 
much beyond the statistical uncertainty of the CA re- 
sults. Correspondingly, we make sure that the biases 
due to a finite time step and a finite number of walk- 
ers are much smaller than the statistical uncertainties of 
the CA results. Furthermore, as discussed in Ref. I23L a 
smoother match to the CA results of the FN energy in 
the /U — » 00 limit is expected using the nodal structure 
given by Slater determinants of plane waves, instead of 
the more accurate*^ (and computationally more demand- 
ing) backflow nodes. 

We simulate N particles in a cubic box with "twist- 
averaged boundary conditions"— (TABC), which have 
been shown to eliminate most of the finite-size effect due 
to the shell structure of the plane-wave determinants. 
For each system considered, simulations are performed 
for 35 points in the irreducible wedge of the first Bril- 
louin zone (BZ) of the simulation box, corresponding to 
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FIG. 1: The size dependence of the VMC energy for the 
Coulomb potential at — 1. Empty symbols: VMC energy 
En; filled symbols; En +Tao—TN ■ The curves show the best- 
fit dependence of the latter, according to Eg. 1251 circles 
and solid line refer to (" = (left scale), triangles and dashed 
line to = 1 (right scale) . The cross is t he C = result ob- 
tained in the thermodynamic limit by Ref. lij. which appears 
fully consistent with the present calculation. The statistical 
errors on the data points are much smaller than the symbol 
sizes. The is 11-8 for C = 0, and 0.4 for C, — l (much poorer 
values would be obtained by just fitting Evmc, i-e., without 
including the kinetic-energy size correction). 



a 1000-point mesh in the whole BZ. 

Both the interparticle potential and the RPA pair 
pseudopotential are computed using an optimized 
breakup^ into a long-range part, to be treated in recip- 
rocal space, and a short-range part, to be treated in real 
space. The short-range part is expanded in locally piece- 
wise quintic Hermite interpolants over 20 knots, and the 
fc-space summation includes 20 shells of reciprocal lat- 
tice vectors. This choice of parameters ensures that, for 
the Coulomb interaction, the potential energy calculated 
for a simulation box containing 64 particles on a simple 
cubic lattice reproduces the exact Madelung constant to 
less than 1 part in 10^. 

All DMC simulations have been done with iV = 54 
for both the paramagnetic and the spin-polarized fiu- 
ids (there is no need of choosing closed-shell determi- 
nants with TABC). Following a common practice )2Sii^ 
the residual size effect has been estimated assuming that 
it is the same for DMC and Variational Monte Carlo 
(VMC), 41 which is somewhat less accurate but much 
cheaper. Systems with up to 246 particles were simu- 
lated with the VMC algorithm, and the size dependence 
of the computed energies was modeled as 



£'00 — E, 



N 



N 



-/3/N, 



(25) 



where T^o and Tjv are the kinetic energy in the thermody- 
namic limit and in the A^— particle system (with TABC), 
respectively, and E^o and /3 are fitting parameters. The 
value, less than 2 on average, is at worst about 10 for 
rs = l and /i = 4, at C = 0. FigureHshows the size extrap- 
olation procedure for the Coulomb potential at = 1. 
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Since the dependence on spin polarization of the optimal 
value of P is very weak (see Fig. QJ, a systematic study 
of the finite-size effect was carried out only for C = 0: for 
given /i and r^, the same value of /3, determined from 
the VMC energies of the paramagnetic fluid at several 
system sizes, was then used to estimate the finite-size 
correction to the DMC energy for both C = and C = 1- 
For the unpolarized gas, we found that the discrepancies 
on the correlation energy with the coupled-cluster data 
of Refs. andlll are of the order of 5-8%. 



IV. ANALYTIC REPRESENTATION OF THE 
CORRELATION ENERGY 

We construct an analytical representation of the cor- 
relation energy as 



to the DMC data. For the parallel-spin g"{0,rs,C) ^tnd 
for the on-top ^(O, rs,C) an exchange-like C dependence 
was assumed, starting from the values at C = l and C = 0, 
respectively. The on-top 17(0, r^,^ = 0) was taken from 
Ref. while g"{0, rs, C = 1) was obtained as a best fit 
to our DMC data. The parameters ai{rs,() of Eq. H26|l 
are then equal to 

ai = 466^3 + 6^^5, 

as = blC3, 

ai = blC2+Utec, 

as = b^ec, 

where ec(^s: C) is the parametrization of the CA correla- 
tion energy as given by Perdew and Wang,"^^ and 



MO 



(26) 



where the function Q is given by Eq. (|22|l . the parameters 
ai{rs, C) ensure the correct large-/i behavior of Eqs. H23|) - 
(|^ . and bo{rs) is fixed by a best fit to our DMC data. 
Some more free parameters which are adjusted to fit the 
DMC data are also contained in the coefficients ai{rs, C), 
whose definition requires a detailed explanation. For the 
limits of Eqs. (|23() - (|24|l we use, for any spin-polarization 
C, the approximation 



ec(»'s,C) - 



3(l-C')gc(0,r„C-0) 
2,g(0,r„C = 0) 9c4(r„C) 



'27r 11^ 



9c5(rs,C) 




C2 
C3 
C4 
C5 



3(l-C^)ffc(0,r„C^0) 
Br? 



2,g(0,r„C = 0) 



9c4(r„C) 

64rf 
9c5(rs,C) 



2tt r5 



40V2^r3 



(30) 



The functions 69, g" , D2 and are finally obtained from 
a best fit to the DMC data and read 



'(0,r„C = l) = 



D2{rs) 



6o(r-s) = 0.784949 r, (31) 

2'/' 1 - 0.02267r, 

5a2r2 (l + 0.4319rs + 0.04r2) ^ ' 

g-0.547r3 

5 (-0.388rs + 0.676r2) (33) 



-0.31rs 



(-4.95r, + rf) . (34) 



The function (f)g, is defined by Eq. H14|l : D2{rs) and I?3(?'s) 
mimic the effect of the ti correlation on the and 
large-/x coefficients, and are obtained by a best fit 



Notice that, by our construction, Eq. H32|) satisfies the ex- 
act high-density limit.*® Our data and the fitting function 
of Eq. (|26|l are shown in Fig. |21 The small discrepancy at 
large /i, particularly visible for — 1 and = 1 on the 
scales of the figure, is due to the condition that our fitting 
function recovers in the Coulomb limit the Perdew- Wang 
parametrization^ of the CA correlation energy, and it is 
consistent with our FN results being an upper bound 
to the data obtained^*' by CA using a nominally exact 
method. 



V. PAIR-DISTRIBUTION FUNCTIONS AND 
ALTERNATIVE SEPARATION OF EXCHANGE 
AND CORRELATION 

From our DMC runs we also extracted, in the usual 
way,~ the pair-distribution functions (/'"^(r, Tg, C, /x). A 
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FIG. 2; Our DMC data for the correlation energy (•) of the 
electron gas with long-range interaction erf (/j,r)/r are com- 
pared with the fitting function (lines) of Eq. 12611 for the 
unpolarized case (upper panel) and the fully polarized case 
(lower panel). The statistical errors on the DMC data are 
comparable with the symbol size. 



FIG. 3: A sample of our DMC pair-distribution functions. 
Upper panel: for rs = 2 and C — 0, g^^ is shown for fi = 
0.25, 0.5, 1, 2, and for the Coulomb gas (/i = oo). Lower panel: 



for rs = 10 and — 1, g is shown for fi - 
and for the Coulomb gas (/i = cxd). 



: 0.05, 0.1, 0.2, 0.4, 



sample of our results is shown in Fig. |31 These func- 
tions are of interest in the framework of the approach 
of Refs. I7''8'l3'l4* 15lllll7l While a local- (or local-spin- 
) density approximation for both exchange and correla- 
tion has been, and to a large extent still is, the most 
popular approach to Kohn-Sham calculations (possibly 
with GGA improvements), there is a growing interest^ 
in optimized-effective-potential schemes, where the ex- 
change is treated exactly and the construction of ap- 
proximations only concerns the correlation energy. The 
latter is naturally defined as whatever exceeds the exact- 
exchange energy, obtained from a single Slater determi- 
nant of Kohn-Sham orbitals. But once a multidetermi- 
nantal, partially correlated wavefunction 'i'^ [Eq. |SJ)] is 
introduced, as in the modified schemes we are concerned 
with here, an alternative, more efficient choice, may be 
to construct approximations only for that portion of the 
correlation energy which is not already taken into ac- 
count by \E'''. In other words, one may prefer to define^S 
"exchange" and "correlation" energy functionals in the 
following way: 

KmdK,^i] = {^^\Vee-Vl^\^n-E'^[n] (35) 
^^c,md["T''^J = Kci^hni]- E';^^^^[n^,ni], (36) 



and then apply, e.g., the LSD approximation only to the 
"correlation" energy functional of Eq. 136|) :^®.. 



■E^c,md["T''^i] = y^^Wcc, md(?-s(r),C(r),M)- (37) 
Here 



the mixed term Alr-sr(^sj Ci m) is equal to 



(38) 



ALR-SR(rs,C,M) /Wdr5LR(r,r„C,M)^^^^^, 
^ Jo ^ 

(39) 

and g];^ is given by minus the pair-distribution func- 
tion of the noninteracting gas. Using the results of Ref.ls^ 
it is easy to show that, for large /z, the mixed term 
Alr-sr behaves as 



Atr_ 



LR-SR 



3 5c(0,r3,C) 



g(0,r„ 0(272-1) 



+ 0{n-') (40) 
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for C 7^ 1; E^nd as 
Alr-sr 



95^'(0,r„C = l) 



/i-+oo 64 rf/i^ 

3g"(0,r„C=l)(3-V2) 



20\/27rr3 



0(^-«)(41) 



for C = li with the same notations of Eas. H23l24() . 

In this Section we present an accurate parametriza- 
tion of Alr-sr- Exploiting our DMC pair-distribution 
functions g^^{r, Vg, Ci A*): we solved Eq. (|39|l by numerical 
integration, and parametrized our results as 



iLR-SR 



(42) 



where the fmictions 5i{rs, C) with i — 3... 6 guarantee the 
correct large- /i behavior of Eqs. (|40|I - H41(I : 



63 = ^dlc^ + dlc^ 

Si = AdlC2+dlCi 
S5 = dlCs 
Se = dlC2. 

Here C2{rs,C) ^nd C4{rs,C) are those of Eqs. (|30|l : 
7(0,r„C = 0)(2\/2-l) 



(43) 
(44) 
(45) 
(46) 



C3 = -(1-n 

^5 = - 



20Fr3 
3c5(r,,C)(3- V2) 



20V2^r3 



(47) 



(7(0,rs,C = 0) and C5(rs,C) are defined in Sec. II VI The 
remaining parameters S2{rs) and do{rs,^) are fitted to 
our DMC data and read 



52{rs) = 0.073867 
do{r,X) = (0.70605 + 0.12927 C^) 

VI. CONCLUSIONS 



(48) 
(49) 



We have presented a comprehensive numerical and an- 
alytic study of the ground-state energy of a homogeneous 
electron gas with modified, long-range-only electron- 
electron interaction erf (fir) /r, as a function of the cutoff 
parameter ^, of the electronic density, and of spin polar- 
ization. The final outcome of this work is the publica- 
tion of a reliable local-spin-density functional which fits 
the results of our quantum Monte Carlo simulations and 
automatically incorporates exact limits. Such a func- 
tional (Sec. IIV|I . or its variant implying the use of an 
additional term also obtained in this work (Sec.0), are 
the key ingredient for some recently proposed "multide- 
terminantal" versions of the density functional theory, 
where quantum chemistry and approximate exchange- 
correlation functionals are combined to optimally de- 
scribe both long- and short-range electron correlations. 



A fortran subroutine that evaluates our LSD 
exchange-correlation functional and the corre- 
sponding potentials is available upon request to 
gori@lct.jussieu.fr, or can be downloaded at 
http : //www. let . jussieu.fr/DFT/gori/elegas .html, 
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APPENDIX A: DETAILS OF EQ. llT^ 

We start from the RPA equations^ for the spin- 
polarized electron gas, and we simply replace the 
Coulomb interaction 1 /r with the long-range interaction 
erf (/ir)/r. We then repeat the analysis done in Refs. Issi 
and 38 for the coulombic gas and find that, in the 
limit, the correlation energy is given by 



12 f°° dy 



„ y / du< a Rciu) e 
I" ./o a Jo 



y 



where 



1 












2 







with 



R{u) = — [1 — uarctan(u "'")] 



(Al) 



(A2) 



(A3) 



zi = (l + C)'/^ Z2 = (l-C)'/^ and a = (4/97r)i/3. 
Equation (|A1I) already shows that the correlation energy 
becomes a function of ix^/rl in the — » limit. 

To prove the small- a:; behavior of Eq. (|20|1 . take a value 
of yuyfrl = a <C 1. In this case all the contribution to the 
integral of Eq. (|A1|) comes from small y, since, as soon as 
?/ 3> a, the integrand goes to zero exponentially fast, as 
a function of a, when a ^ . We thus integrate over y 
Eq. (|A1|) between and a value q\ such that a <C gi <C 1. 
Since y <C 1, the integral reduces to 



/ dyye 



duR({u), (A4) 
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012345678 



X 

FIG. 4: The numerical evaluation of Eq. ijAl^ . as a function of 
X — /iy^/(/)2(C)) a-nd multiplied by [(j)2{C)]~'' ■ If the scaling of 
Eq. 1119^ were exact, all the values corresponding to different 
would lie on the solid curve. The value = 0.86 shown in 
the figure corresponds to the maximum deviation from the 
scaling of Eq. <AH . 



which gives, to leading orders in a when a —^ 0, Eqs. H19(l - 



(|20|l . The large x behavior of Eq. (|21|l follows by consid- 
ering the n ^ oo limit of Eq. I|A1|I . which reduces to the 
standard coulombic case studied in Ref. IHE 

The C dependence of Eq. 1)19(1 is exact in the a: ^ 
limit of Eq. I|A4|) . For larger x, we evaluated Eq. IjAlll nu- 
merically, and in Fig. 0] we report our results multiplied 
by [(f>2{C)]~^: as a function of X = H^/r^/(t'2{C)- if the scal- 
ing of Eq. (|I9|I were exact, all the values corresponding 
to different ( would lie on the solid curve, corresponding 
to C = and 1. The value ( = 0.86 reported in the figure 
corresponds to the maximum deviation from the scaling 
of Eq. (|I9|I . which is thus rather small. The function 
Q{x) of Eq. (1^ has been obtained by fitting the RPA 
data of the solid curve. On the scale of Fig. 01the fitting 
error is invisible. 

To conclude the discussion, we expect that the correla- 
tion energy lies on the curve of Fig. 0] when iiTs <SC 1 
(high-density or really long-range-only interaction on the 
scale Ts). This means that at a given rs, the "exact" 
e];^ lies on the curve of Fig. ^ for values of n such that 
P'\P^s ^ '^1 \f^s^ that is, only the densities < 0.1 would 
be affected by the small deviations from the scaling in C 
of Eq. H19|l . which appear at a; > 3. 



^ W. Kohn, Rev. Mod. Phys. 71, 1253 (1999). 

^ A.E. Mattsson, Science 298, 759 (2002). 

^ C. Fiolhais, F. Nogueira, and M. Marques (eds.), A Primer 

in Density Functional Theory (Springer- Verlag, Berlin, 

2003). 

* In this paper we deal with the separation of the long-range 
and the short-range part of the electron-electron interac- 
tion within DFT, a strategy which can efficiently describe 
van der Waals interactionsii8*lS' For other, recent, promising 
attempts to take into account dispersion forces within DFT 
see, e.g., M. Dion, H. Rydberg, E. Schroder, D. Langreth, 
and B. Lundqvist, Phys. Rev. Lett. 92, 246401 (2004); J. 
Dobson and J. Wang, B. Dinte, K. McLennan, and H. Le, 
Int. J. Quantum Chem. 101, 579 (2005); E. Johnson and 
A. Becke, J. Chem. Phys. 123 024101 (2005); A. Becke and 
E. Johnson, ihid. 122, 154104 (2005); 123, 154101 (2005). 

^ J. Ireta, J. Neugebauer, and M. Scheffler, On the accuracy 

of DFT for describing hydrogen bonds: dependence on the 

bond directionality, preprint. 
® H. Stoll and A. Savin, in Density Functional Method in 

Physics, edited by R. M. Dreizler and J. da Providencia 

(Plenum, Amsterdam, 1985); 

A. Savin, in Recent Developments and Applications of Mod- 
ern Density Functional Theory, edited by J.M. Seminario 
(Elsevier, Amsterdam, 1996). 

* T. Leininger, H. Stoll, H.-J. Werner, and A. Savin, Chem. 
Phys. Lett. 275, 151 (1997); R. PoUet, A. Savin, T. 
Leininger, and H. StoU, J. Chem. Phys. 116, 1250 (2002). 

^ W. Kohn, Y. Meir, and D. E. Makarov, Phys. Rev. Lett. 

80, 4153 (1998). 
" H. Ikura, T. Tsuneda, T. Yanai, and K. Hirao, J. Chem. 

Phys. 115, 3540 (2001); M. Kamiya, T. Tsuneda, and K. 

Hirao, ibid. 117, 6010 (2002); T.T.Y. Tawada, S. Yanagi- 

sawa, T. Yanai, and K. Hirao, ibid. 120, 8425 (2004). 



J. Heyd, G.E. Scuseria, and M. Ernzerhof, J. Chem. Phys. 
118, 8207 (2003). 
^2 R. Baer and D. Neuhauser, Phys. Rev. Lett. 94, 043002 
(2005). 

J. Toulouse, F. Colonna, and A. Savin, J. Chem. Phys. 
122, 014110 (2005). 

R. PoUet, F. Colonna, T. Leininger, H. Stoll, H.-J. Werner, 
and A. Savin, Int. J. Quantum Chem. 91, 84 (2003). 
J. Toulouse, F. Colonna, and A. Savin, Phys. Rev. A 70, 
062505 (2004). 

J. G. Angyan, I.C. Gerber, A. Savin, and J. Toulouse, 
Phys. Rev. A 72, 012510 (2005); I. Gerber and J. G. 
Angyan, Chem. Phys. Lett. 415, 100 (2005). 
1^ E. GoU, H.-J. Werner and H. Stoll, Phys. Chem. Chem. 
Phys. 7 , 3917 (2005). 

S. Yamanaka, K. Kusakabe, K. Nakata, T. Takada, and K. 
Yamaguchi , physics/0508121 

K. Burke, J. P. Perdew, and M. Ernzerhof, J. Chem. Phys. 
109, 3760 (1998); 
^° J.P. Perdew, A. Savin, and K. Burke, Phys. Rev. A, 51, 
4531 (1995). 

2^ M. Levy, Proc. Natl. Acad. Sci. U.S.A. 76, 6062 (1979). 
J. Toulouse and A. Savin, J. Mol. Struct. (Theochem), to 
appear. 

L. Zecca, P. Gori-Giorgi, S. Moroni, and G. B. Bachelet, 
Phys. Rev. B 70, 205127 (2004). 

A. Savin and H.-J. Flad, Int. J. Quantum Chem. 56, 327 
(1995); J. Toulouse, A. Savin, and H.-J. Flad, tbid. 100, 
1047 (2004). 

U. von Barth and L. Hedin, J. Phys. C 5 1629 (1972). 
R. O. Jones and O. Gunnarsson, Rev. Mod. Phys. 61, 689 
(1989). 

2^ A. Nagy, Phys. Rep. 298, 1 (1998). 

^® QMC provides the most accurate numerical values for sev- 



9 



eral physical properties of the HEG. See, e.g., G. Giuliani 
and G. Vignale, Quantum Theory of the Electron Liquid 
(Cambridge, 2005); D. M. Ceperley, in The electron liquid 
paradigm in condensed matter physics, G. F. Giuliani and 
G. Vignale eds. (lOS, Amsterdam, 2004). 
J. Toulouse, P. Gori-Giorgi, and A. Savin, Theor. Chem. 
Acc. 114, 305 (2005). 

D.M. Ceperley and B.J. Alder, Phys. Rev. Lett. 45, 566 
(1980). 

S.H. Vosko, L. Wilk, and M. Nusair, Can. J. Phys. 58, 
1200 (1980). 

J.P. Perdew and A. Zunger, Phys. Rev. B 23, 5048 (1981). 
J.P. Perdew and Y. Wang, Phys. Rev. B 45, 13244 (1992). 
see, e.g., D. Pines and P. Nozieres, Theory of Quantum 
Liquids (Benjamin, New York, 1966). 

Y. Wang and J.P. Perdew, Phys. Rev. B 44, 13298 (1991). 
J.P. Perdew and Y. Wang, Phys. Rev. B 46, 12947 (1992); 
56, 7018(E) (1997). 

P. Gori-Giorgi and J.P. Perdew, Phys. Rev. B 66, 165118 
(2002). 

Y. Wang and J.P. Perdew, Phys. Rev. B 43, 89 11 (1991). 
P. Gori-Giorgi and A. Savin, |cond-mat/051122ri to appear 
in Phys. Rev. A. 

P. Gori-Giorgi and J.P. Perdew, Phys. Rev. B 64, 155102 



(2001). 

For a recent review on fixed-node diffusion Monte Carlo 

and further references, see M. Foulkes, L. Mitas, R. Needs, 

and G. Rajagopal, Rev. Mod. Phys. 73, 33 (2001). 

D. Ceperley, Phys. Rev. B 18, 3126 (1978). 

Y. Kwon, D. M. Ceperley, and R. M. Martin, Phys. Rev. 

B 58, 6800 (1998). 

C. Lin, F. H. Zong, and D. M. Ceperley, Phys. Rev. E 64, 
016702 (2001). 

V. NatoU and D. M. Ceperley, J. Comput. Phys. 117, 171 
(1995). 

V. A. Rassolov, J. A. Pople, and M. A. Ratner Phys. Rev. 
B 62, 2232 (2000). 

see, e.g., G. Ortiz and P. Ballone, Phys. Rev. B 50, 1391 
(1994); 56, 9970(E) (1997). 

see, e.g., S. Kiimmel and J.P. Perdew, Phys. Rev. B 68, 
035103 (2003); W. Yang and Q. Wu, Phys. Rev. Lett. 89, 
143002 (2002); R. J. Magyar, A. Fleszar, and E. K. U. 
Gross, Phys. Rev. B 69, 045111 (2004); M. Griining, O. 
V. Gritsenko, and E. J. Baerends, J. Chem. Phys. 118, 
7183 (2003). For a critical review, see also E.J. Baerends 
and O.V. Gritsenko, J. Chem. Phys. 123, 062202 (2005). 



